Predicting Financial Crashes Using Discrete Scale Invariance 



Anders Johansen 1 , Didier Sornette 1 ' 2,3 and Olivier Ledoit 4 
1 Institute of Geophysics and Planetary Physics 
University of California, Los Angeles, California 90095 

2 Department of Earth and Space Science 
University of California, Los Angeles, California 90095 
3 Laboratoire de Physique de la Matiere Condensee 
CNRS UMR6622 and Universite de Nice-Sophia Antipolis 
B.P 71, Pare Valrose, 06108 Nice Cedex 2, France 
4 Anderson Graduate School of Management 
University of California, Box 90095-1481, 110 Westwood Plaza 
Los Angeles CA 90095-1481 

February 1, 2008 



Abstract 

We present a synthesis of all the available empirical evidence in the light of recent theoretical 
developments for the existence of characteristic log-periodic signatures of growing bubbles in a 
variety of markets including 8 unrelated crashes from 1929 to 1998 on stock markets as diverse as 
the US, Hong-Kong or the Russian market and on currencies. To our knowledge, no major financial 
crash preceded by an extended bubble has occurred in the past 2 decades without exhibiting such 
log-periodic signatures. 



1 Introduction 



Two recent papers [ Johansen and Sornette, 1999a , Johansen et al, 2000 ] have presented increasing 
evidence on the Oct. 1929, Oct. 1987, Hong-Kong Oct. 1987 crashes, on the Aug. 1998 global market 
events and on the 1985 Forex event on the US dollar, for the following hypothesis [Sornette et al, 1996] 
describing stock market crashes : we propose that they are caused by a slow build-up of long-range time 
correlations reflecting those between traders leading to a collapse of the stock market in one critical 
instant. This build-up manifest itself as an over- all power law acceleration in the price decorated by 
"log-periodic" precursors. Here, "log-periodicity" refers to a sequence of oscillations with progressively 
shorter cycles of a period decaying according to a geometrical series. In addition, extensive statistical 
tests have been performed [ pohansen et al., 200C , Johansen, 1997] ] to show that the reported "log- 
periodic" structures essentially never occurred in ~ 10 5 years of synthetic trading following a "classical" 
time-series model, the GARCH(1,1) model with student-t statistics, often used as a benchmark in 
academic circles as well as by practitioners. Thus, the null hypothesis that log-periodicity could result 
simply from random fluctuations is strongly rejected. 

From a theoretical view point, we have proposed a rational expectation model of bubbles and 
crashes, which has two main components : (1) We assume that a crash may be caused by local self- 
reinforcing imitation processes between noise traders that can be quantified by the theory of critical 
phenomena developed in the Physical Sciences ; (2) We allow for a remuneration of the risk of a crash 
by a higher rate of growth of the bubble, which reflects that the crash is not a certain deterministic 
outcome of the bubble and, as a consequence, it remains rational for traders to remain invested 
provided they are suitably compensated. 

Since our initial proposition [fSornette et al, 1996 ], several works have extended the empirical inves- 
tigation [Feigenbaum and Freund, 199€, Feigenbaum and Freund, 1998, Gluzman and Yukalov, 1998, 
Vandcwalle et al, 1998a, Vandewalle et al, 1998b, pohansen and Sornette, 1999b , Drozdz et al., 1999| , 
Johansen and Sornette, 1999c | as well as the theoretical framework of critical phenomena 
| Sornette and Johansen, 1997, [Sornette and Johansen, 1998 ], However, a recent article of opinion 
[ Laloux et al, 1999 1 and a short assessment note [ Dinski, 199S [ have raised concerns about the re- 
ality of the log-periodic signatures in bubbles preceding crashes and the relevance of the concept of 
critical phenomena to financial crashes in general. 

The purpose of the present paper is primarily to present a complete up-to-date synthesis of the 
available evidence containing all previously reported cases together with four new cases which include 
the correction of the US dollar against the Canadian dollar and the Japanese Yen starting in Aug. 
1998, as well as the bubble on the Russian market and its ensuing collapse in 1997-98. We are able 
to show a remarkable universality of the results for all events, with approximately the same value of 
the fundamental scaling ratio A characterising the log-periodic signatures. We also address briefly the 
concern raised by Ilinski | Ilinski, 1999 ] about the validity of our model | Johansen and Sornette, 1999a , 
Johansen et al., 200C]. Considering the robustness of the results presented here, our analysis opens 
the road to financial crash forecasting (see for example footnote [12] of | Stauffer and Sornette, 1998 1 
and jVandewalle et al, 1998bfl ). 

A similar analysis of log-periodicity for "anti-bubbles" have been used to issue such a forecast 
in Dec. 1998 for the recovery of the Nikkei in 1999 [Johansen and Sornette, 1999c ]. At the time of 
writing (May 1999), the forecast, performed at a time when the Nikkei was at its lowest, seems to 
have correctly captured the change of regime and the overall upward trend since the beginning of this 
year. 
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2 Theoretical framework 



In order to put the empirical analysis in perspective, it is useful to recapitulate the main features of 
the model [Johansen and Sornette, 1999a, Johansen et al., 2000| 1 . This will also allows us to stress its 



robustness with respect to other choices in the model construction and address the criticism made by 



Ilinski HHinski, 1999 | 1 



2.1 Ingredients 

• Our key assumption is that a crash may be caused by local self-reinforcing imi- 
tation between traders [Shiller, 1990, Shiller, 1995, Censi et al, 1996, Cutler et al, 1990| , 



Keim and Madhavan, 19*95] , Nelson, 1994 ] . This self-reinforcing imitation process leads to the 
blossoming of a bubble. If the tendency for traders to "imitate" their "friends" increases up 
to a certain point called the "critical" point, many traders may place the same order (sell) at 
the same time, thus causing a crash. The interplay between the progressive strengthening of 
imitation and the ubiquity of noise requires a probabilistic description : A crash is not a certain 
outcome of the bubble but can be characterised by its hazard rate h{t), i.e., the probability per 
unit time that the crash will happen in the next instant provided it has not happened yet. 

• Since the crash is not a certain deterministic outcome of the bubble, it remains rational for 
traders to remain invested provided they are compensated by a higher rate of growth of the 
bubble for taking the risk of a crash, because there is a finite probability of "landing smoothly" , 
i.e., of attaining the end of the bubble without crash. In this model, the ability to predict the 
critical date is perfectly consistent with the behaviour of the rational agents: They all know this 
date, the crash may happen anyway, and they are unable to make any abnormal risk-adjusted 
profits by using this information. 

Our model distinguishes between the end of the bubble and the time of the crash : The rational 
expectation constraint has the specific implication that the date of the crash must have some degree 
of randomness. The theoretical death of the bubble is not the time of the crash because the crash 
could happen at any time before, even though this is not very likely. The death of the bubble is only 
the most probable time of the crash. 

The model does not impose any constraint on the amplitude of the crash. If we assume it propor- 
tional to the current price level, then the natural variable is the logarithm of the price. If instead, we 
assume that the crash amplitude is a finite fraction of the gain observed during the bubble, then the 



natural variable is the price itself [ Johansen and Sornette, 1999a ] . We are aware that the standard 



economical proxy is the logarithm of the price and not the price itself, since only relative variations 
should play a role. However, as we shall see below, different price dynamics gives both possibilities. 

In the stylised framework of a purely speculative asset that pays no dividends and in which we 
ignore the interest rate, risk aversion, information asymmetry, and the market-clearing condition, 
rational expectations are simply equivalent to the familiar martingale hypothesis: 

for all t' > t E t [p(t')} =p{t) , (1) 

where p(t) denotes the price of the asset at time t and Et[-] denotes the expectation conditional on 
information revealed up to time t. In the sequel, we will relax the hypothesis of no risk aversion and 
show that the predictions of the model are very robust. This is important as risk aversion is often 
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considered as a key ingredient controlling the crash dynamics and the time-dependence of the risk 
perception by market participants also controls in part the pre-crash dynamics. 

We model the occurrence of a crash as a jump process j whose value is zero before the crash and 
one afterwards. The random nature of the crash occurrence is modeled by the cumulative distribution 
function Q(t) of the time of the crash, the probability density function q(t) = dQ/dt and the hazard 
rate h{t) = q(t)/[l — Q(t)]. The hazard rate is the probability per unit of time that the crash will 
happen in the next instant provided it has not happened yet. 

Assume for simplicity that, during a crash, the price drops by a fixed percentage k G (0,1), say 
between 20 and 30% of the price increase above a reference value p±. Then, the dynamics of the asset 
price before the crash are given by: 



dp = fj,(t) p(t) dt - K[p(t) - Pl ]dj , (2) 

where j denotes a jump process whose value is zero before the crash and one afterwards. Taking Ej [dp] 
and using the fact that Ejd/] = h(t)dt leads to 

H{t)p{t) = K\p{t) - Pl ]h{t) . (3) 

In words, if the crash hazard rate h(t) increases, the return \x increases to compensate the traders for 
the increasing risk. Plugging (§) into (|2[), we obtain a ordinary differential equation. For p(t) —p(to) < 
p(to) — Pi-, its solution is 



p(t) ^ p(t ) + K[p(t ) - px] [ h(t')dt' before the crash. (4) 

Jtn 



The integral f to h(t')dt' is the cumulative probability of a crash until time t. We have found that this 
regime, where the price itself is the relevant observable, applies to the relatively short time scales of 
approximately two to three years prior to the crash studied here. 

If instead the price drops by a fixed percentage k € (0, 1) of the price, the dynamics of the asset 
price before the crash is given by 

dp = n(t)p(t) dt — np(t)dj . (5) 

We then get 

E t [dp] = fi(t)p(t)dt - Kp(t)h(t)dt = , (6) 

which yields : 

fi(t) = Kh(t) . (7) 
and the corresponding equation for the price is : 

p(t) 1 rt 



log 



p(*o) 



k / h(t')dt' before the crash. (8) 
Jt 



This gives the logarithm of the price as the relevant observable. It has successfully been ap- 
plied to the 1929 and 1987 Wall Street crashes up to about 7.5 years prior to the crash 
[Sornette and Johansen, 1997, Johansen et al, 200C]. 



The higher the probability of a crash, the faster the price must increase (conditional on having no 
crash) in order to satisfy the martingale condition. Intuitively, investors must be compensated by a 
higher return in order to be induced to hold an asset that might crash. This is the only effect that 
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we wish to capture in this part of the model. This effect is fairly standard and it was pointed out 
earlier in a closely related model of bubbles and crashes under rational expectations by Blanchard 
( fBlanchard, 197SH , top of p.389). It may go against the naive preconception that price is adversely 
affected by the probability of the crash, but our result is the only one consistent with rational expec- 
tations. Notice that price is driven by the hazard rate of crash h(t). 

Ilinski [Ilinski, 1999] raises the concern that the martingale condition (||) leads to a model which 



"assumes a zero return as the best prediction for the market." He continues : "No need to say that 
this is not what one expects from a perfect model of market bubble! Buying shares, traders expect 
the price to rise and it is reflected (or caused) by their prediction model. They support the bubble 
and the bubble support them!". 

In other words, Ilinski (1999) criticises a key economic element of our model : Market rationality. 
We have captured this by assuming that the market level is expected to stay constant as written in 
equation ([!]). Ilinski (1999) claims that this equation (|l]) is wrong because the market level does not 
stay constant in a bubble: It rises, almost by definition. 

This misunderstanding addresses a rather subtle point of the model and stems from the difference 
between two different types of returns : 

1. The unconditional return is indeed zero as seen from (jlj) and reflects the fair game condition. 

2. The conditional return, conditioned upon no crash occurring between time t and time t' , is non- 
zero and is given by equations (||) and (|7|), respectively. If the crash hazard rate is increasing 
with time, the conditional return will be accelerating precisely because the crash becomes more 
probable and the investors need to be remunerated for their higher risk. 

Thus, the expectation which remains constant in equation (Q) takes into account the probability 
that the market may crash. Therefore, conditionally on staying in the bubble (no crash yet), the 
market must rationally rise to compensate buyers for having taken the risk that the market could have 
crashed. 

The market price reflects the equilibrium between the greed of buyers who hope the bubble will 
inflate and the fear of sellers that it may crash. A bubble that goes up is just one that could have 
crashed but did not. Our model is well summarised by borrowing the words of another economist : 
"(...) the higher probability of a crash leads to an acceleration of [the market price] while the bubble 
lasts." Interestingly, this citation is culled from the very same article by Blanchard [ Blanchard, 1979| 



that Ilinski [Ilinski, 1999] cites as an alternative model more realistic than ours. We see that this is 
in fact more of an endorsement than an alternative. 

To go into details, Blanchard's model [ Blanchard, 1979| is slightly more general because it incor- 



porates risk aversion through : 

uE t \p(t')} = p(t) , (9) 

where v £ (0, 1] is an appropriate discount factor. We have been aware of this since the beginning, 
and the only reason why we did not take it into account is that it obviously makes no difference in 
our substantive predictions (log-periodic oscillations and power law acceleration), as long as v remains 
bounded away from zero and infinity. 

Another way to incorporate risk aversion into our model is to say that the probability of a crash in 
the next instant is perceived by traders as being K times bigger than it objectively is. This amounts 
to multiplying our hazard rate h(t) by K, and once again this makes no substantive difference as long 
as K is bounded away from zero and infinity. 
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The point here is that v and K both represent general aversion of fixed magnitude against a risk. 
Risk aversion is a central feature of economic theory, and it is generally thought to be stable within 
a reasonable range, associated with slow-moving secular trends such as changes in education, social 
structures and technology. 



Ilinski | Ilinski, 1999 1 rightfully points out that risk perceptions are constantly changing in the 
course of real-life bubbles, but wrongfully claims that our model violates this intuition. In our model, 
risk perceptions do oscillate dramatically throughout the bubble, even though subjective aversion to 
risk remains stable, simply because it is the objective degree of risk that the bubble may burst that goes 
through wild swings. 

For these reasons, the criticisms put forth by Ilinski, far from making a dent in our economic 
model, serve instead to show that it is robust, flexible and intuitive. 

2.2 Crash hazard rate and critical imitation model 

The crash hazard rate h{t) quantifies the probability that a large group of agents place sell orders 
simultaneously and create enough of an imbalance in the order book for market makers to be unable to 
absorb the other side without lowering prices substantially. Most of the time, market agents disagree 
with one another, and submit roughly as many buy orders as sell orders (these are all the times when 
a crash does not happen). The key question is to determine by what mechanism did they suddenly 
manage to organise a coordinated sell-off? 



We have proposed the following answer [ Johansen et al, 2000 ]: All the traders in the world are 



organised into a network (of family, friends, colleagues, etc.) and they influence each other locally 
through this network : For instance, an active trader is constantly on the phone exchanging information 
and opinions with a set of selected colleagues. In addition, there are indirect interactions mediated for 
instance by different parts of the media. Specifically, if I am directly connected with k other traders, 
then there are only two forces that influence my opinion: (a) The opinions of these k people and of 
the global information network and (b) an idiosyncratic signal that I alone generate. Our working 
hypothesis here is that agents tend to imitate the opinions of their connections. The force (a) will tend 
to create order, while force (b) will tend to create disorder. The main story here is a fight between order 
and disorder. As far as asset prices are concerned, a crash happens when order wins (everybody has 
the same opinion: Selling), and normal times are when disorder wins (buyers and sellers disagree with 
each other and roughly balance each other out). This mechanism does not require an overarching 
coordination mechanism since macro-level coordination can arise from micro-level imitation and it 
relies on a realistic model of how agents form opinions by constant interactions. 

Many models of interaction and imitation between traders have been developed. To make a long 
story short, the upshot is that the fight between order and disorder often leads to a regime where order 
may win. When this occurs, the bubble ends. Many models (but not all) undergo this transition in a 
"critical" manner [Goldenfeld, 1992, Dubrulle et al, 1997] : The sensitivity of the market reaction to 



news or external influences accelerate on the approach to this transition in a specific way characterized 
by a power law divergence at the critical time t c of the form F(t) = (t c — t)~ z , where z is called a 
critical exponent. This form amounts to the property that 

dlnF 



d1n(t c - t) 



(10) 



is a constant, namely that the behavior of the observable F becomes self-similar close to t c with respect 
to dilation of the distance t c — t to the critical point at t c . The symmetry of self-similarity in the 
present context refers to the fact that the relative variations dlnF = dF/F of the observable with 
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respect to relative variations dln(t c — t) = d(t c — t)/{t c — t) of the time-to-crash are independent of 
time t, as expressed by the constancy of the exponent z. 

Accordingly, the crash hazard rate follows a similar dependence, namely 



h(t) = B(t c - t) 



-b 



(11) 



where B is a positive constant and t c is the critical point or the theoretical date of the bubble death. 
The exponent b must lie between and 1 for an important economic reason : otherwise, the price 
would go to infinity when approaching t c (if the bubble has not crashed yet). 

We stress that t c is not the time of the crash because the crash could happen at any time before t c , 
even though this is not very likely. t c is the most probable time of the crash. There exists a residual 
probability 



of attaining the critical date without crash. This residual probability is crucial for the coherence of 
the story, because otherwise the whole model would unravel because rational agents would anticipate 
the crash. 



where = 1 — b G (0, 1) and p c is the price at the critical time t c (conditioned on no crash having 
been triggered). We see that the price before the crash also follows a power law with a finite upper 
bound p c . The slope of the price, which is the expected return per unit of time, becomes unbounded 
as we approach the critical date t c . This is to compensate for an unbounded hazard rate approaching 



In the contemporary information environment, one could argue that the media can play a more 
efficient role in the creation of an opinion than the network of "friends" and will thus modify if not 
destroy the critical nature of the crash resulting from the increasing cooperativity between the network 
of "friends". There are two ways to address this issue. First, consider the situation of a typical trader. 
Since all her competitors have access to the same globally shared information, she does not feel that 
this gives her an edge, but only serves to keep her abreast the market moves. The informations that 
she consider really valuable are those that she can share and exchange with a selected and narrow 
circle of trusted colleagues. This argument brings us back to the model of an imitation network with 
"local" connections. A second argument is that the media do nothing but reflect a kind of consensus 
resulting precisely from the collective action of all the actors on the market. It is thus a kind of 
revelator of the interactions between the traders. In addition, this revelation of the general consensus 
and of the global market sentiment reinforces the effective interaction between the traders. In the 
language of the Statistical Physics of critical points, the media can thus be compared to an effective 
mean-field created endogenously and may thus effectively strengthen the existence of a critical point 
as described here. 

2.3 Log-periodicity 

The power law dependence ( [TT|) of the hazard rate is the hallmark of self-similarity of the market 
across scales as the critical time t c is approached : at the critical point, an ocean of traders who are 
mostly bearish may have within it several islands of traders who are mostly bullish, each of which in 




(12) 



Plugging equation ([Tl]) into Equation (||) gives the following law for price: 

kB 

p(t) rs p c 7r{tc — i) before the crash. 



(13) 
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turns surrounds lakes of bearish traders with islets of bullish traders; the progression continues all the 
way down to the smallest possible scale: a single trader [Wilson, 1979]. Intuitively speaking, critical 



self-similarity is why local imitation cascades through the scales into global coordination. This critical 
state occurs when local influences propagate over long distances and the average state of the system 
becomes exquisitely sensitive to a small perturbation, i.e., different parts of the system becomes highly 
correlated. 

As we said, scale invariance of a system near its critical point must be represented by power law 



dependences of the observables Dubrulle et al, 1997 ] . Formally, these power laws are the solution of 



the renormalisation group equations [Wilson, 1979, Goldenfeld, 1992] which describe the correlations 



between agents at many scales. It turns out that the most general mathematical solutions of these 
renormalisation group equations are power laws with complex exponents which as a consequence 
exhibits log-periodic corrections to scaling |Sornette, 199Sfl . 

A straightforward mechanism for these complex exponents to appear is to define the dynamics on 
a hierarchical structure displaying a discrete scale invariance. Schematically, we can think of the stock 
market made of actors which differs in size by many orders of magnitudes ranging from individuals to 
gigantic professional investors, such as pension funds. Furthermore, structures at even higher levels, 
such as currency influence spheres (US$, Euro, YEN ...), exist and with the current globalisation and 
de-regulation of the market one may argue that structures on the largest possible scale, i.e., the world 
economy, are beginning to form. This means that the structure of the financial markets have features 
which resembles that of hierarchical systems with "traders" on all levels of the market. Of course, 
this does not imply that any strict hierarchical structure of the stock market exists, but there are 
numerous examples of qualitatively hierarchical structures in society. Models [ lohansen et al., 2000| , 



Sornette and Johansen, 1998(1 °f imitative interactions on hierarchical structures recover the power 
law behavior (^) and predict that the critical exponents b and (3 may be complex numbers! 

Recently, it has been realized that discrete scale invariance and its associated complex expo- 
nents may appear "spontaneously" in non-hierarchical systems, i.e., without the need for a pre- 
existing hierarchy (see [ [Sornette, 1998 [ for a review). There are many dynamical mechanisms 



that can operate to produce these complex exponents, such as non-local geometry, frozen het- 
erogeneity Saleur and Sornette, 1996|], cascades of instabilities, intermittent amplifications, (see 



[ Johansen and Sornette, 1999d] for such a mechanism in the context of "anti-bubbles" proposed re- 



cently for the Japanese stock market), and so on. 

As a consequence, we are led to generalize ( |TT|) and give the first order expansion of the general 
solution for the hazard rate 

h(t) « B(t c - t)~ b + C(t c - t)- b cos[uilog(t c - t) + if)]. (14) 

The crash hazard rate now displays log-periodic oscillations. This can easily seen to be the general- 
ization of equation ( |l0|) by taking the exponent z to be complex with a non-zero imaginary part, since 
the real part of (t c — t)~ z+tuJ is (t c — t)~ z cos (u ln(t c — £))• The evolution of the price before the crash 
and before the critical date is then given by: 

pit) « A 1 + Bi(i c - tf + d(t c - tf cos (wi log(t c - t) + epi) (15) 

where <j)\ is another phase constant. The key feature is that accelerating oscillations appear in the 
price of the asset before the critical date. The local maxima of the function are separated by time 
intervals that tend to zero at the critical date, and do so in geometric progression such that the ratio 
of consecutive time intervals is a constant 
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A = e 2 ^ 1 . (16) 

This is very useful from an empirical point of view because such oscillations are much more strik- 
ingly visible in actual data than a simple power law : a fit can "lock in" on the oscillations which 
contain information about the critical date t c . 

We note that these log-periodic structure bear some similarity with patterns classified em- 
pirically by chartists and others using methods of technical analysis, such as "Elliott waves" 



[ Frost and Prechter, 1998 1 and "log-spirals" | Erman, 1999 ] . For instance, a logarithmic spiral in the 
plane obeys the equation r = tq e ad in polar coordinate. The intersections with the x axis occur 
for 6 = nir (recall that x = rcosO) and the corresponding coordinates are x n = (— l) n e awn . The 
intersections with the y axis occur for 9 = tt/2 + nir (recall that y = rsinO) and the corresponding 
coordinates are y n = (—1) e an ^ 2 e ann . Both series indeed form discrete geometrical series with a 
scaling ratio A = e a7T similar to (|l6|). We do not need, however, to stress the difference between a 



pattern recognition approach | Erman, 1999 ] and our economically motivated approach incorporating 



a rational model of imitative agents as well as well-established tools from statistical physics. 

3 Empirical results 

3.1 Are large crashes special? 

Crashes are extreme events. There are two possibilities to describe them : 

1. The distribution of returns is stationary and the extreme events can be extrapolated as lying 
in its far tail. Within this point of view, recent works in finance and insurance have recently 
investigated the relevance of the body of theory known as Extreme Value Theory to extreme 
events and crashes [ Embrechts et al., 1997 , Embrechts et al., 1998 , Embrechts et al., 1998j ], 



2. Crashes cannot be accounted for by an extrapolation of the distribution of smaller events to the 
regime of extremes and belong intrinsically to another regime, another distribution, and are thus 
outliers. 

To test which one of these two alternatives provides the most accurate description, a statistical 
analysis of market fluctuations [ Johansen and Sornette, 1998| ] has provided significant indications that 



large crashes may be outliers. Indeed, it was established that the distribution of draw downs of the Dow 
Jones Average daily closing is well described by an exponential distribution with a decay constant of 
about 2%. This exponential distribution holds only for draw downs smaller than about 15%. In other 
words, this means that all draw downs of amplitudes of up to approximately 15% are well-described 
by the same exponential distribution with characteristic scale 2%. This characteristic decay constant 
means that the probability to observe a draw down larger than 2% is about 37%. Following hypothesis 
1 and extrapolating this description to, e.g., the 3 largest crashes on Wall Street in this century (1914, 
1929 and 1987) yields a recurrence time of about 50 centuries for each single crash. In reality, the 
three crashes occurred in less than one century. This suggests that hypothesis 2 is preferable. In the 
following, we use the generic term "crash" to refer to the significant drop of a market occurring over 
a few days to a few weeks which follows an extended period of bullish behaviour and signifies the end 
of the bubble. 

As an additional null-hypothesis, 10.000 synthetic data sets, each covering a time-span close to a 
century hence adding up to about 10 6 years, have been generated using a GARCH(1,1) model estimated 
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from the true index with a t-student distribution with four degrees of freedom [ Bollerslev et al, 1992 ], 
This model is often used in academic and practitioner circles as a reasonable description of the market 
and includes both the non-stationarity of volatilities and the fat tail nature of the price returns. Our 



analysis | Johansen et al., 2000| ] showed that only two data sets had 3 draw downs above 22% and 



none had 4. However, 3 of these 6 "crashes" were preceded by a draw up of comparable size of the 
draw down and hence showed an abnormal behaviour not found for real crashes. This means that in 
approximately one million years of "Garch-trading" , with a reset every century, never did 3 crashes 
occur in a single century. 

This suggest that different mechanisms may be responsible for large crashes and that hypothesis 
2 is the correct description of crashes. This also raises the issue whether pre-cursory patterns exists 
decorating the speculative bubble preceding the crash. The point is simply that while the GARCH(1,1) 
model does a reasonable job of reproducing fluctuations in "normal trading", it cannot capture the 
fluctuations connected with large crashes and hence another type of model is necessary for this special 
behaviour. Of course, these simulations do not prove that our model is the correct one, only that one 
of the standard models of the "industry" (which makes a reasonable null hypothesis) is utterly unable 
to account for the stylized facts associated to large financial crashes. A better model is thus called for 
and this is exactly the goal of the present paper to provide for such a framework. 

3.2 Log-periodic analysis of large crashes 

The numerical procedure of fitting equation (|l5|) , as well as equations (|l^) , (|2(]) and ( pl| ) is a minimi- 
sation of the variance 

1 N 

V*r = -Y,{yi-f{ti)f (17) 

i=l 

between the ./V data points y« and the fit function / (t) using the down-hill simplex algorithm 
[ Press et al, Tggl . Tests using the maximum likelihood method with the student-t distribution show 



that our results are not significantly sensitive to the estimation procedure [ Lin, 1998 ]. In order to 
reduce the complexity of the fit, all linear variables A, B, C are determined by the nonlinear variables 
by requiring that a solution has zero derivative of the variance as a function of the linear variables. 
This means that equation (15), is an effective 4-parameter fit with f},<f>,t c ,u. Further details on our 



numerical procedure are given in [ Johansen et al., 200C| , Johansen, 1997]. 



A remaining question concerns what data interval prior to the crash to fit. The procedure used 
was the following. The last point used for all crashes was the highest value of the price before the 
crash, the first point used was the lowest value of the price when the bubble started. For all the cases 



discussed here, except the Russian anti-bubble analysed in section |3.5f| , this procedure always gave a 
convincing result and it was never necessary to change the interval fitted. 

Figures [l] to || show the behaviour of the market index prior to the 4 stock market crashes of Oct. 
1929 (Wall Street), Oct. 1987 (Wall Street), Oct. 1997 (Hong-Kong) and Aug. 1998 (Wall Street) as 
well as the collapse of the US$ against the DEM and CHF in 1985 and against the Canadian dollar 
and the YEN in 1998. A fit with equation (|l^) is shown as a continuous line for each event. Table |j] 
lists the parameters of the fit to the data. Note the small fluctuations in the value of the scaling ratio 
2.2 < A < 2.8 for all data sets except the CHF. This agreement cannot be accidental and constitutes 
one of the key test of our framework. 

In order to qualify further the significance of the log-periodic oscillations in a non-parametric way 
as well as providing an independent test on the value of the frequency of the log-periodic oscillations 

*Here "lowest" and "highest" is of course interchanged 
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crash 


tc 


^max 




drop 


P 




A 


A 


B 


C 


Var 


1929 (WS) 


30.22 


29.65 


29.87 


47% 


0.45 


7.9 


2.2 


571 


-267 


14.3 


56 


1985 (DEM) 


85.20 


85.15 


85.30 


14% 


0.28 


6.0 


2.8 


3.88 


1.16 


-0.77 


0.0028 


1985 (CHF) 


85.19 


85.18 


85.30 


15% 


0.36 


5.2 


3.4 


3.10 


-0.86 


-0.055 


0.0012 


1987 (WS) 


87.74 


87.65 


87.80 


30% 


0.33 


7.4 


2.3 


411 


-165 


12.2 


36 


1997 (H-K) 


97.74 


97.60 


97.82 


46% 


0.34 


7.5 


2.3 


20077 


-8241 


-397 


190360 


1998 (WS) 


98.72 


98.55 


98.67 


19% 


0.60 


6.4 


2.7 


1321 


-402 


19.7 


375 


1998 (YEN) 


98.78 


98.61 


98.77 


21% 


0.19 


7.2 


2.4 


207 


-84.5 


2.78 


17 


1998 (CANS) 


98.66 


98.66 


98.71 


5.1% 


0.26 


8.2 


2.2 


1.62 


-0.23 


-0.011 


0.00024 



Table 1: i c is the critical time predicted from the fit of the market index to the equation (|l5|). The 
other parameters of the fit are also shown. The error Var is the variance between the data and the 
fit, equation (|i~7|), and has units price 2 . Each fit is performed up to the time t max at which the market 
index achieved its highest maximum before the crash. t m i n is the time of the lowest point of the 
market disregarding smaller "plateaus". The percentage drop is calculated from the total loss from 

tmax to tmin- 



and hence the prefered scaling ratio A, we have eliminated the leading trend from the price data by 
the transformation 

p{t)-A 1 -B 1 {t c -tf 

P{t) ^ d(t c -ty • (18) 



This transformation should produce a pure cosfwi log(t c — t) + <f>i] if equation (|15| ) was a perfect 
description. In figure |j| we show the residual for the bubble prior to the 1987 crash. We see a very 
convincing periodic behaviour as a function of log (~f~^^J ■ I n figure |(| we see the corresponding data 
as a function of the true time to the crash t c — t. Looking at figures || and |(| especially two things 
strikes the eye. First, we see that in the latter figure the distance between consecutive extrema of the 
data is clearly accelerating, in the former we clearly see that a transformation to the logarithm of the 
time to the crash exhibits an over-all periodic trend, hence demonstrating a log-periodic behaviour. 
Second, since the detrended data oscillates nicely around zero, the two figures clearly show that the 



power law in equation (15) is an excellent approximation to the over-all rise in the data. That the 
residual oscillates beyond ±1 only indicates that the cosine used in equation ( |T5| ) is not an optimal 
choice for a periodic function in the present context. However, it has the advantage of being familiar 
to the majority as well as straightforward to implement numerically. 

The significance of this oscillating trend has been further estimated by performing a spectral 
analysis, designed to test for the statistical significance of the log-periodic oscillations. In the goal, we 



use a so-called Lomb periodogram [Press et al, 1994] rather than the more standard power spectrum 



since the data is not equidistantly sampled in time. The Lomb periodogram is a local fit of a cosine 
(with a phase) using some user chosen range of frequencies. In figure [7|, we see a peak around / ~ 1.1 
for all 8 cases shown in figures |l| to f| corresponding to u\ = 2nf ~ 7 in perfect agreement with the 
results listed in table |l[ Since the "noise" spectrum is unknown and very likely different for each 
crash, we cannot estimate the confidence interval of the peak in the usual manner [press et al, 1994| 
and compare the results for the different crashes. Therefore, only the relative level of the peak for 
each separate periodogram can be taken as a measure of the significance of the oscillations and the 
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periodograms have hence been normalised. We note that if the noise were to have a white Gaussian 
distribution, the confidence given by the periodograms increases exponentially on the value of the 



abscissa and would be well above 99.99% for all cases shown | Press et al, 1994 1. Note also, that the 
strength of the oscillations is ~ 5 — 10% of the leading power law behaviour for all 8 cases signifying 
that they cannot be neglected. 

To summarize thus far, our spectral analysis demonstrates that the observed log-periodic oscilla- 
tions have a very strong power spectrum, much above noise level. It would be very difficult and much 
less parsimonious to account for these structures by another model. 

3.3 Are log-periodic signatures statistically significant? 

In the case of the 1929 and 1987 crashes on Wall Street, log-periodic oscillations in the index could 



be identified as long as ~ 7.5 years prior to the crashes Sornette and Johansen, 1997 1. In order to 
investigate the significance of these results, we picked at random fifty 400-week intervals in the period 
1910 to 1996 of the Dow Jones average and launched the same fitting procedure as on the time period 
prior to the 1929 and 1987 crashes on Wall Street. The results were encouraging. Of 11 fits with a 
quality (specifically, the variance of the fit of equation ( pO] ) with r = t c — t to the data) comparable with 
that of the 2 crashes, only six data sets produced values for (3, u\, u>2 and T\, which were in the same 
range. However, all of these fits belonged to the periods prior to the crashes of 1929, 1962 and 1987. 
The existence of a "crash" in 1962 was before these results unknown to us and the identification of this 
crash naturally strengthens the case. We refer the reader to [Johansen et al, 200C] for a presentation 
of the best fit obtained for this "crash" . 

We also generated 1000 synthetic data sets of length 400 weeks using the same GARCH(1,1) model 



used in section |3.l| . Each of these 1000 data sets was analysed in the same manner as the real crashes. 
In 66 cases, the best minima had values similar to the real crashes, but all of these fits, except two, 
did not resemble the true stock market index prior to the 1929 and 1987 crashes on Wall Street, the 
reason primarily being that they only contained one or two oscillations. However, two fits looked 



rather much like that of the 1929 and 1987 crashes [Johansen et al, 2000, Johansen, 1997]. 

This result correspond approximately to the usual 95 % confidence interval for one event. In 
contrast, we have here provide 5 examples of log-periodic signatures before large stock market crashes. 
We are aware that this is not proof of log-periodic signatures in bubbly markets, this only a thorough 
data analysis can only establish within reasonable doubt, but it provides for a very good estimation 
that these signatures are not easily generated accidentally. 

3.4 Log-periodicity in bearish markets? 

Stock market jargon divide the stock markets trends into either "bullish" or "bearish". If log-periodic 
signals, apparently, are present in "bullish" markets, the obvious question to ask is whether this is 
also the case in "bearish" markets. 

The most recent example of a genuine long-term depression comes from Japan, where the Nikkei 
has decreased by more than 60 % in the 9 years following the all-time high of 31 Dec. 1989. In figure 
we see (the logarithm of) the Nikkei from 31 Dec. 1989 until 31 Dec. 1998. The 3 fits are equations 




and ( pi|) respectively [ Johansen and Sornette, 1999c | : 



log (p(t)) » Ai + B^ 13 + dr p cos [u! log (r) + fa] (19) 
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log (p(t)) 



A 2 + 



1 + 



2/1 



B 2 + C 2 cos 



U) 2 

u\ log r + — log 



1 + 



+ 



(20) 



log (p(t)) « A3 + 



4£i 



-B3 + C3 cos 



+ 03 



(21) 



where r = t — t c . Equation (|2l]) predicts the transition from the log- frequency u\ close to t c to 



u>i + u>2 for T\ < t < T2 and to the log-frequency u)\ + 0J2 + for T2 < r. The equations (|20j ) (resp. 
(|2l|)) extend the renormalisation group approach to the second (resp. third) order of perturbation 



Johansen and Sornette, 1999c , Sornette and Johansen, 1997 1 . 



The interval used for equation ( |19| ) is until mid-1992 and for equation (pQ ) until mid-1995. The 
parameter values produced by the different fits with equations (^) and (|20|) agree remarkably well 
and the values for the exponent (3 also agree well as shown in table ^. The fit with (|2l"|) was, due to 
the large number of free variables, performed differently and the parameter values for t c , /3, and u\ 
determined from the fit with (p0|) was kept fixed and only T\, T2, u>2, <^3 and <ps where allowed to 



adjust freely. What lends credibility to the fit with equation (21) is that despite it complex form, we 
get values for the two cross-over time scales T\, T2 which correspond very nicely to what is expected 
from the theory: T\ has moved down to ~ 4.4 years in agreement with the time interval used for 
equation (20) and T2 is ~ 7.8 years, which is compatible with the 9 year interval of the fit. This does 
not mean that fitting is not very degenerate, it is, but the ranking of T% and T2 is always the same 
and the values given do not deviate much from the ones in the caption of figure ^ i.e., by ± one year. 

The value uj\ ~ 4.9 correspond to A ~ 3.6, which is not in the range found for the bubbles, see 
table 0. An additional difference between the anti-bubble in the Nikkei and these other cases is the 
strength of the oscillations compared to the leading behaviour. For the Nikkei, it is ~ 20%, i.e., 2 — 4 
times as large as the amplitude obtained for the stock market and Forex bubbles. 



3.5 Log-periodic analysis of the Russian stock market 



Intrigued by the claims of K. Ilinski [ Ilinski, 1999H of a genuine stock market bubble, followed by 



a crash, without log-periodic signatures on the Russian stock market, we have analysed 4 Russian 
stock market indices from their start up till present. The 4 indices was The Russian Trading System 
Interfax Index (IRTS), The Agence Skate Press Moscow Times Index (ASPMT), The Agence Skate 
Press General Index (ASPGEN) and The Credit Suisse First Boston Russia Index (ROSI). The ROSI 
is generally considered the best of the 4 and we will put emphasis on the results obtained with that 
index. As we shall return to later, the Russian stock market is highly volatile, which means that great 
care must be taken in maintaining a representative stock market index. This is the primary reason 
for using 4 indices in the analysis. 

In figure H we see the ROSI fitted with equation ( [To]) in the interval [96.21 : 97.61]. The interval 
is chosen by identifying the start of the bubble and the end represented by the date of the highest 
value of the index before the crash similarly to the 7 major market crashes discussed previously. For 
all 4 indices, the same start- and end-day could be identified ±1 day. 

In figures and H we see the detrended data using the transformation (|l8|). The conclusions 
made in relation to the corresponding figures for the 1987 Wall street crash (figures || and ^) are again 
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valid supporting both an over-all power law rise as well as log-periodic signatures for the Russian crash 
of 1997. 

As can be seen from table ^, the non-dimensional parameters j3, uj\ and A as well as the predicted 
time of the crash t c for the fit to the different indices agree very well except for the exponent (3 obtained 
from the ASPGEN Index. In fact, the value obtained for the prefered scaling ratio A is fluctuating by 
no more than 5% for the 4 fits showing numerical stability. 



Bubble 


tc 


tmax 


tmin 


drop 


P 




A 


A 


B 


C 


Var 


ASPMT 


97.61 


97.61 


97.67 


17% 


0.37 


7.5 


2.3 


1280 


-1025 


59.5 


907 


IRTS 


97.61 


97.61 


97.67 


17% 


0.39 


7.6 


2.3 


633 


-483 


38.8 


310 


ROSI 


97.61 


97.61 


97.67 


20% 


0.40 


7.7 


2.3 


4254 


-3166 


246 


12437 


ASPGEN 


97.62 


97.60 


97.67 


8.9% 


0.25 


8.0 


2.2 


2715 


-2321 


72.1 


1940 


Anti-bubble 


t c 


tmax 


tmin 


drop 


P 




A 


A 


B 


C 


Var 


ROSI 


97.72 


97.77 


98.52 


74% 


0.32 


7.9 


2.2 


4922 


-3449 


472 


59891 


Nikkei eq.( 


m 


89.99 


90.00 


92.63 


63% 


0.47 


4.9 


3.6 


10.7 


-0.54 


-0.11 


0.0029 


Nikkei eq.( 


^cjD 


89.97 


90.00 


95.51 


63% 


0.41 


4.8 


3.7 


10.8 


-0.70 


-0.11 


0.0600 



Table 2: t c is the predicted time of the crash from the fit of the market index to the equation (|l~5|). 
The other parameters of the fits to the preceding bubble are also given. The error Var is the variance 
between the data and the fit, equation ([17]), and has units price 2 except for the Nikkei, where the 
units are [log(price)] 2 . The fit to the bubble is performed up to the time at which the market index 
achieved its highest maximum before the crash. The parameters t c , (3, lo\ and A correspond to the fit 



with equation (15), where t c and t has been interchanged. Here t max and t m i n represent the endpoints 



of the interval fitted. 

The origin of this bubble is well-known. In 1996 large international investors (read US, German 
and Japanese) began to invest heavily in the Russian markets believing that the financial situation 
of Russia had finally stabilised. Nothing was further from the truth |Intriligator, 1998| , |Malki, 1999] 



but the belief and hope in a new investment haven with large returns led to herding and bubble 
development. This means that the same herding, which created the log-periodic bubbles on Wall 
Street (1929, 1987, 1998), Hong Kong (1997) and the Forex (1985, 1998), entered an emerging market 
and brought along the same log-periodic trends characterising the global markets. The fact that the 
consistent values of A obtained for the 4 indices of the Russian market are comparable to that of the 
Wall Street, Hong Kong and Forex crashes supports this interpretation. Furthermore, it supports the 
idea of the stock market as a self-organising complex system of surprising robustness. 

Inspired by this clear evidence of log-periodic oscillations decorating the power law acceleration 
signaling a bubble in the Russian stock market, we extended our analysis and search for possible 
log-periodic signatures in the "anti-bubble" that followed the log-periodic bubble described above. 



As described in section 3.4, the recent decay of the Japanese Nikkei index starting 1. Jan 1990 

we 



until present can be excellently model by a log-periodically decorated power law decay. In figure 10 



show the ROSI index for the "anti-bubble", signaling the collapse of the Russian stock market, fitted 



with equation (15) where t c and t has been interchanged. We use the price and not the logarithm of 
the price to keep the symmetry with the description of the bubble and from the consideration of the 
relatively short time scales involved. 



In table H, the corresponding values for the fit of the anti-bubble with equation (|19|) listed. The 
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agreement between the bubble and the anti-bubble with respect to the values of the non-dimensional 
variables (3, u\ and A as well as t c is surprisingly good. However, the numerical stability of the result 
for the anti-bubble cannot be compared with that of the preceding bubble and depends on the index 
as well as the endpoint of the interval used in the fitting. However, the "symmetry" around t c is rather 
striking considering the nature of the data and quite visible to the naked eye. 

There are (at least) two reasons for the numerical instability of the fit to the anti-bubble. First, 
as an emerging market in decay the Russian stock market is too volatile for a smooth function as a 
cosine and we are currently investigating alternative parametrisations [|Johansen and Sornette, 1999d ] . 
Second, the "noise" in the Russian indices were very large in that period, in particular due to the 
"depletion" of stocks. This is clearly illustrated by the heavy rearrangement that the ASP GEN and 
related indices went through following Aug. 17, 1998. 

In figure [l3|, we see the spectral analysis for the bubble and anti-bubble. Whereas the agreement 
between the periodogram and the fit with equation (^) is excellent (within 1%), the deviation for 
the anti-bubble is ~ 5%. Furthermore, the residual for the anti-bubble used in the frequency analysis 
shown in figure |l^ was truncated by a few weeks due to an increasing deviation between data and fit. 

It may seem odd to argue for the log-periodic oscillations while one can forcefully argue that the 
market is largely reflecting the vagaries of the Russian political institutions. For instance, in the 
anti-bubble case, Feb-April 1998 was a revival period for the market characterized by the returning of 
western investors after the post-crash calm-down. This can be followed by studying the dynamics of 
the Russian external reserves. The timing of the return can be argued to be dictated by the risk policies 
of larger investors more than anything else. The next large drop of the Russian index in April 1998 
originated by the decision of Mr. Yeltsin to sack Mr. Chernomyrdin's government, which destabilized 
the political situation and created uncertainty. Further political disturbance was introduced twice by 
the Duma when it rejected Mr. Yeltsin's candidates for the prime minister office and put itself on the 
brink of dissolution. Opposed to this, we argue that one must not mistake a global unstable situation 
for the specific historical action that triggered the instability. Consider a ruler put vertically on a 
table. Being in an unstable position, the stick will fall in some direction and the specific air current 
or slight initial imperfection in the initial condition are of no real importance. What is important 
is the intrinsically unstable initial state of the stick. We argue that a similar situation applies for 
crashes. They occur because the market has reached a state of global instability. Of course, there will 
always be specific events which may be identified as triggers of market motions but they will be the 
revelators rather than the deep sources of the instability. Furthermore, the political events must also 
be considered as revelators of the state of the dynamical system which includes the market. There is, 
in principle, no decoupling between the different events. Specifically, the Russian crash may have been 
triggered by the Asian crises, but it was to a large extent fueled by the collapse of a banking system, 
which in the course of the bubble had created an outstanding debt of $19.2 billion |Malki, 1999 j . 



4 Conclusion 

We have presented a synthesis of the available empirical evidence in the light of recent theoretical 
developments for the existence of characteristic log-periodic signatures of growing bubbles in a variety 
of stock markets as well as currencies. We have here documented 8 unrelated crashes from 1929 to 
1998, on stock markets as diverse as the US, Hong-Kong or the Russian market and on currencies. In 
addition, we have discovered a significant bubble on Wall Street ending in 1962 [ Johansen et a/., 2000H 
as well as "anti-bubbles" on the Nikkei since 1990 and the Gold (after the 1980 bubble maximum) 



Johansen and Sornette, 1999c ] . Quite unexpectedly, we have shown that the Russian bubble crashing 
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in Aug. 1997 had close to identical power law and log-periodic behaviour to the bubbles observed on 
Wall Street, the Hong-Kong stock market and on currencies. To our knowledge, no major financial 
crash preceded by an extended bubble has occurred in the past 2 decades without exhibit a log- 
periodic signature. In this context, note that the novel analysis of the Russian index presented here 
was motivated by Ilinski's claim of a crash without log-periodic signature, which we have shown to be 
incorrect. 

All these results, taken together with the remarkable robustness and consistency of the estimation 
of the exponent (3 as well as the more important statistics the scaling ratio A, make the case for 
power law acceleration and log-periodicity very strong. In our opinion, one can no more ignore theese 
very specific and strong signatures which is characteristic of developing bubbles and this calls for 
further investigations to unravel in more depths the underlying economical, financial and behavioural 
mechanisms. 

These different cases, together with the 1962 "slow event" as well as the "anti-bubbles", show that 
the log-periodic critical theory applies both to bubbles ending in a sudden crash as well as to bubbles 
landing smoothly. This is in fact a strong prediction of our rational model of imitative behaviour. 

What we have attempted here is not to explain why crashes happens or bubbles exists, but to 
quantify the process taking place during extended bubbles that very often lead to the rapid regime- 
switching a crash represent. We have offered evidence that bubbles and anti-bubbles have log-periodic 
and power law characteristics and hence provides for a quantification of extended "moods" on the 
markets. We have furthermore, provided a model which contains the observed signatures. Whether, 
as proposed, discrete scale invariance or some other mechanisms [Sornette, 1995] are responsible for 
the prominent log-periodic signatures observed, only a more detailed microscopic model can answer. 
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Figure 1: The Dow Jones index prior to the October 1929 crash on Wall Street. The fit is equation 
@ with A 1 w 571, Si ps -267, C x ~ 14.3, /? « 0.45, t c w 30.22, <^i » 1.0, wi « 7.9. 
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Figure 2: The S&P 500 US index prior to the October 1987 crash on Wall Street and the US $ 
against DEM and CHF prior to the collapse mid-85. The fit to the S&P 500 is equation fll5| ) with 
Ai rj 412, Bi « -165, d w 12.2, (3 w 0.33, t c w 87.74, fa ps 2.0, wi » 7.4. The fits to the DM and 
CHF currencies against the US dollar gives A\ w 3.88, B x w -1.2, Ci « 0.08, /3 ps 0.28, t c ps 85.20, 
01 ps -1.2, uji ps 6.0 and Ai » 3.1, Bi ps -0.86, d ps 0.05, /? ps 0.36, t e w 85.19, 0i ps -0.59, 
wi ps 5.2, respectively. 
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Figure 3: The Hang Seng index prior to the October 1997 crash on the Hong-Kong Stock Exchange 
and the S&P 500 stock market index prior to the recent crash on Wall Street in August 1998. The 
fit to the Hang Seng index is equation © with A x ps 20077, B x ps -8241, |B X ps -397, (3 ps 0.34, 
t c ps 97.74, 0i ps 0.78, wi ps 7.5. The fit to the S&P 500 index is equation fl]]| with A 1 ps 1321, 
5i « -402, Ci ps 19.7, (3 ps 0.60, t c ps 98.72, 0i ps 0.75, ^ « 6.4. 
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Figure 4: The CAN$ and YEN currencies against the US dollar prior to the drop starting in Aug. 
1998. The fit with equation (15) to the two exchange rates gives Ai ps 1.62, B\ ps —0.22, C\ ps —0.011, 
(3 ps 0.26, t c ps 98.66, 0i ps -0.79, o>i ps 8.2 and A-y ps 207, Bi ps -85, d ps 2.8, /3 ps 0.19, t c ps 98.78, 
4>\ ps —1.4, wi ps 7.2, respectively. 
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Figure 5: The residual as denned by the transformation (|l^) as a function of log f^pM for the 1987 
crash. 
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Figure 6: The residual as defined by the transformation ([i8|) as a function of t c — t for the 1987 crash. 
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Figure 7: The Lomb periodogram for the 1929, 1987 and 1998 crashes and Wall Street, the 1997 crash 
on the Hong Kong Stock Exchange, the 1985 US $ currency crash in 1985 against the DM and CHF 
and in 1998 against the YEN and the 5.1% correction against the CAN$. For each periodogram, the 
significance of the peak should be estimated against the noise level. 
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Figure 8: Natural logarithm of the Nikkei stock market index after the start of the decline 1. Jan 
1990 until 31 Dec. 1998. The lines are equation ( |T9| ) (dotted line) fitted over an interval of ~ 2.6 years, 
equation (|^) (continuous line) over as 5.5 years and equation ( |2l|) (dashed line) over 9 years. The 
parameter values of the first fit of the Nikkei are A\ ~ 10.7, B\ ~ — 0.54, Ci ~ — 0.11, /? ~ 0.47, t c ~ 
89.99,01 ~ — 0.86,u;i ~ 4.9 for equation (19). The parameter values of the second fit of the Nikkei 
are A 2 « 10.8, B 2 ^ -.70, C 2 « -0.11, /3 fa 0.41, t c w 89.97, ^ 2 « 0.14, wi w 4.8, T x w 9.5,^2 ~ 4.9 
for equation (|20|). The third fit uses the entire time interval and is performed by adjusting only T\, 



T 2 , oj 2 and oj^ , while j3, t c and wi are fixed at the values obtained from the previous fit. The values 
obtained for these four parameters are T\ ~ 4.3 years, T 2 ~ 7.8 years, oo 2 ~ —3.1 and T 2 ~ 23. Note 
that the values obtained for the two time scales T\ and T 2 confirms their ranking. This last fit predicts 
a change of regime and that the Nikkei should increase in 1999. 
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Figure 10: Symmetric "bubble" and "anti-bubble" : in addition to the ascending part of the ROSI 
Index which is reproduced from figure || with the same fit, we show the deflating part fitted with 
equation [l|. The parameter values of the fit are A ± « 4922, B 1 -3449, C\ « 472, j3 « 0.32, t c « 
97.72,01 « 1.4, wi w 7.9. 
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Figure 11: The residual as defined by the transformation (|lq ) as a function of log f ^ — J for the 
Russian data shown in figure ||. 
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Figure 12: The residual as defined by the transformation ( J18| ) as a function of t c — t for the Russian 
data shown in figure |9| 
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Figure 13: The Lomb periodogram for the bubble and anti-bubble shown in figure 10. The significance 
of the peak should be estimated against the noise level. 
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